A Dynamical Solution to the Problem of a Small Cosmological Constant and 

Late-time Cosmic Acceleration 



C. Armendariz-Picon, 1 V. Mukhanov 1 and Paul J. Steinhardt 2 
1 Ludwig Maximilians Universitat, Sektion Physik, Miinchen, Germany 
2 Department of Physics, Princeton University, Princeton, N J 0854-0 



o ■ 

O ' 

o : 

(N ■ 
< 



> ■ 

-xf : 

m . 
■ 

o : 
o . 
o ■ 

^ ' 

Or 

6 ■ 
a: 



Increasing evidence suggests that most of the energy den- 
sity of the universe consists of a dark energy component with 
negative pressure, a "cosmological constant" that causes the 
cosmic expansion to accelerate. In this paper, we address the 
puzzle of why this component comes to dominate the uni- 
verse only recently rather than at some much earlier epoch. 
We present a class of theories based on an evolving scalar field 
where the explanation is based entirely on internal dynamical 
properties of the solutions. In the theories we consider, the 
dynamics causes the scalar field to lock automatically into 
a negative pressure state at the onset of matter-domination 
such that the present epoch is the earliest possible time, con- 
sistent with nucleosynthesis restrictions, when it can start to 
dominate. 



Introduction. Observations of large scale structure, 
searches for Type la supernovae, and measurements of 
the cosmic microwave background anisotropy all suggest 
that the universe is undergoing cosmic acceleration and 
is dominated by a dark energy component with negative 
pressure.D The dark energy may consist of a cosmologi- 
cal constant (vacuum density) or quintessence, □ such as a 
scalar field with negative pressure. In either case, a key 
challenge is the "cosmic coincidence" problem: Why is 
it that the vacuum density or scalar field dominates the 
universe only recently? Until now, either cosmic initial 
conditions or model parameters (or both) had to be tuned 
to explain the low density of the dark energy component. 

In this paper, we explore a new class of scalar field 
models with novel dynamical properties that avoid the 
fine-tuning problem altogether. A feature of these 
models is that the negative pressure results from the 
non-linear kinetic energy of the scalar field, which we 
call, for brevity, /c-field or /c-essence. (This consider- 
ation is inspired by earlier studies of /c-inflation, ki- 
netic energy driven inflationtrla) . As we will show, 
for a broad class of theories, there exist attractor solu- 
tions which determine the equation-of-state of A:-essence 
during different epochs depending on the equation-of- 
state of the background. Effectively, the scalar field 
changes its speed of evolution in dynamic response to 
changes in the background equation-of-state. During the 
radiation-dominated epoch, fc-essence is led to be sub- 
dominant and to mimic the equation-of-state of radi- 
ation. Hence, the ratio of fc-essence to radiation den- 
sity remains fixed. When the universe enters the dust- 



dominated epoch, though, /c-essence is unable to mimic 
the dust-like equation-of-state for dynamical reasons. In- 
stead, the energy decreases rapidly by several orders of 
magnitude and freezes at a fixed value. After a pe- 
riod (typically corresponding roughly to the current age 
of the universe), the field overtakes the matter density 
and drives the universe into cosmic acceleration. Ulti- 
mately, the fc-essence equation-of-state slowly relaxes to 
an asymptotic value between and -1. (The reader may 
wish to sneak a peek at Fig. 3 which illustrates the be- 
havior in a specific numerical example.) 

The scenario bears some likeness to the quintessence 
"tracker models" discussed by Zlatev et aZ.,t§0 but is 
in fact very different from them. For a certain class 
of tracker potentials, the quintessence scalar field con- 
verges to an attractor solution in which the energy den- 
sity in the quintessence field mimics the equation-of-state 
of the background (matter or radiation) energy density. 
However, for tracker potentials, it does not make a dif- 
ference what the equation-of-state of the background is. 
Only after the scalar field passes a certain critical value 
does the quintessence develop a negative pressure, and, 
then, the energy density becomes fixed. The weakness in 
this model is that the energy density for which the pres- 
sure becomes negative is set by an adjustable parameter, 
which has to be fine tuned to explain why cosmic accel- 
eration is happening at the present epoch in the history 
of the universe. 

The distinctive feature of /c-essence models is that 
tracking of the background energy density can only occur 
in the radiation epoch. At the matter-radiation equal- 
ity, a sharp transition of fc-essence from positive to neg- 
ative pressure is automatically triggered by dynamics. 
The ^-essence cannot dominate before matter-radiation 
equality because it is exactly tracking the radiation back- 
ground. It also cannot dominate immediately after dust- 
domination because its energy density necessarily drops 
several orders of magnitude at the transition to dust- 
domination. However, since its energy density decreases 
more slowly than the matter density as the universe ex- 
pands, fc-essence must dominate not too long thereafter, 
at roughly the current epoch. The resolution of the cos- 
mic coincidence problem boils down to the fact that we 
live at the "right time" after matter-radiation equality. 

As noted above, the remarkable behavior comes at the 
cost of introducing a non-linear kinetic energy density 
functional of the scalar field and adjusting it to obtain 
the desired attractor behavior. This kind of action may 
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describe a fundamental scalar field or be a low-energy 
effective action. For example, in string and supergrav- 
ity theories, non-linear kinetic terms appear generically 
in the effective action describing moduli and massless 
degrees of freedom (superpartners) due to higher_prder 
gravitational corrections to the Einstein action.Em The 
attractor behavior of our models relies on certain broad 
conditions on the form of these terms. Our initial exam- 
ples are admittedly contrived for the purposes of numer- 
ical illustration. A systematic study ..of model-building 
will appear in a forthcoming paperjlj although, hav- 
ing seen here the relatively simple basic principles, the 
reader should be equipped to explore more attractive and 
better-motivated forms. 

Equations. In the theories we consider the Lagrangian 
density for ip is taken to be 



C = --R+^p k (X) + C ri 



(1) 



where R is the Ricci scalar, X = \{Wp) 2 , C m is the 
Lagrangian density for dust and radiation and we use 
units where 8irG/3 = 1. The energy density of the fc-field 
p is pk — (2Xp t x — p)/<y9 2 ; the pressure is pk = p/p 2 ; and 
the speed of sound of fc-essence is c 2 — Pk,x / 'pk ,x jr-frfiere 
the subscript means derivative with respect to XBc 

The attractor behavior can be explained most easily by 
changing variables from X to y = 1/y/X and rewriting 
the fc-field Lagrangian as: 



C k =p k (X)/p 2 = g(y)/tp 2 y. 



(2) 



In this case, the energy density and pressure are pk — 
—g'/p 2 andpfe = g/p 2 y, where prime indicates derivative 
with respect to y. The equation-of-state is 



w k = Pk/pk = -g/yg' 



and the sound speed is 



g'y 



Pk 



g"y 2 



(3) 



(4) 



In order to have a sensible, stable theory, we require 
Pk > and c 2 > 0. These conditions are satisfied if g' < 
and g" > in the region where p' k is positive. There- 
fore, a general, convex, decreasing function g(y), such as 
shown in Fig. 1, satisfies these necessary conditions. Us- 
ing the Friedmann equation: H 2 = p to t = Pfc + Pm , where 
p m is the energy density of matter (radiation and dust), 
and the energy conservation equations, p\ — —3pi(l+Wk) 
for the fc-essence (i = fc) and matter (i = m) components, 
we obtain the following equations of motion 



V= 2 



3 (w k (y) - 1) 



(y) 



r(y)-\l-^- 

Ptot, 



— ) = 3— ( 1 - — ) (w m - w k (y)) , 

yPt0t J Pt0t \ Pt0t , 



(5) 
(6) 



where 

r(y) 
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y (1 + Wk) 



3 is -g'y) 
2V2 V 3 ? : 



(7) 



and dot denotes derivative with respect to TV = In a. 
These are the master equations describing the dynamics 
of fc-essence models. Once some general properties of 
g(y) are specified, the attractor behavior described in 
the introduction follows from these coupled equations. 

Dynamics. We are seeking a tracker solution y(N) 
in which the fc-essence equation-of-state is constant 
and exactly equal to the background equation-of-state, 
Wk(y (N)) = w m , and the ratio pk/ptot is fixed. Generi- 
cally, this requires y(N) be a constant y tr and therefore 
Pk/Ptot = r 2 (y tr ). The last condition can only be sat- 
isfied if r(y tr ) is less than unity. Hence, given a convex 
function such as shown in Fig. 1, we can first identify 
those ranges of y where r(y) is greater than unity or less 
than unity. In ranges where r(y) is less than unity, we can 
seek values of y where Wk in Eq. (||) is equal to w m , the 
equation-of-state of the matter or radiation. The value 
of ytr changes depending on the epoch and w m . These 
are the attractor solutions. In ranges where r(y) exceeds 
unity, there are no attractor solutions. 




FIG. 1. A plot of g(y) vs. y (see Eq. (3) for definition) 
indicating the points discussed in the text. R corresponds to 
the attractor solution during the radiation-dominated epoch; 
S is the de Sitter attractor at the onset of matter-domination; 
and K is the attractor as fc-essence dominates. For our range 
of g(y), there is no dust-like attractor solution at y = yo- 

Although there is no dust attractor, it is quite possi- 
ble for there to be a radiation attractor. The radiation 
attractor corresponds to positive pressure, so it can be 
located only at y < yr>- Hence, we must have g{y) such 
that r(y) is less than unity for some range at y < yo 
which includes some point yn where Wk{yii) = 1/3. Dur- 
ing the radiation-dominated epoch, the ratio of fc-essence 
to the total density remains fixed on this attractor and 
equal to {pk/ Ptot) = r 2 (y R ). 

In Fig. 1, the pressure pk = g/p 2 y is positive above the 
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y-axis, and negative below the y-axis. The dust equation- 
of-state pk = can only be obtained at y = ya where 
g(y) goes through zero. However, this point can be an 
attractor, only if, the second condition, r(yo) < 1, is 
satisfied. If it so happens that r(yn) > 1, then there is 
no dust attractor in the matter-dominated epoch. This 
is precisely what we want for our scenario, and this is 
possible for a broad class of functions g. 

If g possesses a radiation attractor but no dust attrac- 
tor, what happens at dust-radiation equality? To answer 
this question let us study the solutions of the master 
equations Eqs. (|||7|) in two limiting cases, when the en- 
ergy density of fc-essence is either much smaller or much 
greater than the matter energy density. If Pk/pm 
one can neglect the last term in the equation (|5|) and it is 
obvious that y (N) ~ ys, where ys satisfies the equation 
r (ys) — 0, is an approximate solution of the equations 
of motion. The point S satisfies g(ys) = g'(ys)ys, so the 
tangent of g at ys passes through the origin, as shown 
in Fig. 1. Since r oc (1 + itffe), the equation of state 
of fc-essence at ys corresponds to wu (ys) ~ — 1; we call 
this solution the de Sitter attractor and denote it by S 
in Fig. 1. From Fig. 1, it is clear that ys nearly always 
exists for convex decreasing functions g. We stress that 
y (N) ss ys is an approximate solution to the equations- 
of-motion only when matter strongly dominates over fc- 
essence. So, if pk during the radiation dominated epoch 
is significantly less than the radiation density which it 
tracks, which is both typical and required to satisfy nucle- 
osynthesis constraints, then fc-essence proceeds to the de 
Sitter attractor immediately after dust-radiation equal- 
ity.t3 

As the transition to dust-domination occurs, pk first 
drops to a small, fixed value, as can be simply under- 
stood. Suppose that (pk/ptot) R = r 2 (y R ) = a < 1CT 2 
during the radiation dominated epoch, where the bound 
is set by nucleosynthesis constraints. From the equa- 
tion of state, Eq. (||), we have the relation: g(y R ) = 
—g' (y R )y R /3. The condition r(yp) > 1 is required in or- 
der to have no dust attractor solution. Combining these 
relations, we obtain: 

44 < 9 < 10-2 (8) 

9dVd 16 

On the other hand, it is apparent from Fig. 1 that 
—g' R > —g'u, so yji < j/d if a < 1. In particular, 
the tangent at y D falls below g(y R ), so g' D (y R - yn) ~ 
-yog'r, < 9(Vr) = -Vr9r/3- Using this relation, we 
obtain 

M< 3 <2 .l(r 3 and 4 Z <4<7-10- 4 . (9) 
Vd 16 g R 16 

Since p k = -g'/(p 2 and \g' (y s )\ < W (vd)\ we con- 
clude that after radiation domination, when the fc-field 
reaches the vicinity of the S-attractor, the ratio of en- 
ergy densities in fc-essence and dust does not exceed 



(Pk/ptot)n'xg' D /g R ; that is, p k /pdust < a 2 /l6 < 7-10 6 . 
Hence, provided (pk/ptot)R < 10~ 2 at dust-radiation 
equality, the fc-essence field loses energy density on its 
way to the S-attractor down to a value below 7 x 10~ 6 . 

By definition, the S-attractor is one in which w « — 1 
and the energy density is nearly constant. Hence, once 
Pk has reached its small but non-zero value, it freezes. 
In the further evolution of the universe, the matter den- 
sity decreases, but the fc-essence energy density remains 
constant, eventually overtaking the matter density of the 
universe. Note that, as p k approaches p m , the condition 
Pk/ p m < 1 is necessarily violated and a new attractor 
solution is found for the case where fc -essence itself dom- 
inates the background energy density. This attractor is 
denoted K in Fig. 1. 

To prove that that the K-attractor exists, we con- 
sider the master equations, Eqs. in the limit where 
Pk/ptot — > 1. If yK satisfies the equation r (yx) — 1, then 
y (N) ~ yx is an approximate solution of the equations 
of motion. When dust is not a tracker, there always prists 
a unique attractor yx in the interval yn < y < j/s-" To 
prove this, note that, within this interval, the function 
r(y) has a negative derivative. Recall that r(ys) = 
(definition of S-attractor) and r(yo) > 1 (to avoid a 
dust attractor) . Since r (y) is a monotonically decreas- 
ing, continuous function, there exists a unique point yx 
(yn < yK < ys)) where r(y) becomes equal to unity. 
At y > yu the pressure of fc-essence is negative. Hence, 
generically the K attractor, located near yx, describes 
a universe dominated by a negative pressure component 
which induces power-law cosmic acceleration. As accel- 
eration proceeds, pk increasingly dominates and y yx- 

Following along using Fig. 1, the dynamics can be 
summarized as follows: fc-Essence is attracted to y = 
y R during the radiation dominated epoch; at matter- 
domination, the energy density drops sharply as fc- 
essence skips past y = yD, because there is no dust at- 
tractor, and heads towards y w ys- The energy density 
Pk freezes and, after a period, overtakes the matter den- 
sity. As it does so, y relaxes towards yx ■ In this scenario, 
our current universe would be making the transition from 
ys to yx- All this occurs for generic g(y) satisfying broad 
conditions on its first and second derivatives. If the ra- 
tio of pk to the radiation density is near the maximum 
allowed by nucleosynthesis (roughly equipartition initial 
conditions), the scenario predicts that the pk dominates 
by the present epoch. 

Numerical results. We have verified these analytic pre- 
dictions numerically for a wide class of g(y). As as a 
strategy, we look for forms which are roughly linear, 

g (v) « -\g' R yR + 9r(v -vr) + o (fc, - VR ) 2 ) (10) 

in the vicinity of radiation attractor R and parabolic 

9(y)~4 m 2(y-y D )(y-^-)+- (n) 

Vd ~ Vs V Vd) 
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in the region yo < y < ys- One can easily check, that 
the points y&, yn and ys here are, by construction, the 
places where the corresponding attractors are located and 
9ri9'd are the derivatives of g at the appropriate points. 
The results are not sensitive to the precise form of g 
that interpolates between these regimes. The main con- 
straints are that the attractor solution have a small ratio 
of pk/ptot during the radiation epoch and that there is 
no dust attractor. 
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FIG. 2. The ratio of fc-essence energy density, p^, to the 
density in radiation and matter, p m , vs. red shift. At the 
present epoch (dashed line), Qk ~ 0.7. 
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FIG. 3. The fc-essence equation-of-state vs. red 
shift. The three attractors in the radiation-, matter-, and 
fc-essence-dominated epochs are evident. At the present 
epoch, Wk ~ —0.77. 

For illustrative purposes, we have used these principles 
to obtain a sample g{y) and have transformed it to a form 
for p(X) in the action, Eq. (Q): 

Pk(X) = -2.01 + 2VTTX + 0.03(aA:) 3 - {bXf, (12) 

where a = 1CP 5 and b = 10~ 6 . For small values of X, af- 
ter a field redefinition, this Lagrangian density reduces to 
one equivalent to a canonical scalar field with exponential 
potential (curiously, a tracker modeH). The distinctive 
dynamical attractor in our models relies deviations from 



linearity at large X. 

The results of a numerical integration are presented in 
Figs. 2 and 3. We see that fc-essence tracks the radiation 
(wk ~ 1/3) during the radiation-dominated epoch. Then, 
at the onset of matter-domination, u>k starts to change 
and the energy density of fc-essence suddenly drops by 
several orders of magnitude becoming of the order of 10 -7 
times the critical density at red shifts about z ~ 1000 as 
the S-attractor is approached and w — ► — 1. At about red 
shift z ~ 3 — 5, pk becomes non- negligible and Wk starts 
to increase, ultimately reaching Wk — —0.77 at z = 0. 
The ratio of the fc-essence energy density to the critical 
density today is flk ~ 0.74. In the future, Wk in this 
model has to approach the value —0.55, corresponding 
to the K-attractor solution, and the universe will enter 
the period of power law fc-inflation. 

Summary. In this paper, we have presented a scenario 
in which cosmic acceleration occurs late in the history 
of the universe due to an inevitable sequence of events 
caused by attractor dynamics. We view the present work 
as a demonstration of principle; hence, we have empha- 
sized general conditions and an analytic understanding 
of the scenario. The specific example illustrated in this 
paper is admittedly complex, composed to illustrate the 
concept, but we know of no fine-tuning or other require- 
ment that poses a barrier to finding simpler and better- 
motivated forms. By changes of variables and other sim- 
ple techniques, one can quickly enlarge the class of ac- 
tions considered here, which may suggest other types of 
attractive models. 

A prediction of fc-essence models that distinguishes 
them from models based on tracker potentialsQ is that Wk 
is in the process of increasing today from -1 towards its 
asymptotic value at the K attractor, whereas, for track- 
ers, Wk is undergoing a transition from w ps towards 
W = — 1. A consequence is that the effective value of Wk 
for fc-essence - that is, the weighted average of Wk be- 
tween the present and z = 1 - can be significantly lower 
than for the tracker potential case, which is bounded be- 
low by w e ff ps —0.75.0 In the numerical example above, 
the effective w e ff — —0.84, for example. The current 
supernovae data suggest a lower value of Wk more con- 
sistent with fc-essenceJid Of course, the fc-essence range 
for Wk is more difficult to distinguish from a cosmological 
constant [w = — 1)- 

In future work,£3 we discuss model-building using fur- 
ther examples and generalizations. We also explore inter- 
esting variations of the dynamical scenario with different 
kinds of attractors, including some which can lead to dif- 
ferent long-term future outcomes, such as a return to 
a pressureless, unaccelerated expansion in the long-term 
future. 
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